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Abstract 

Bounds for the correlation functions of identical bosons are discussed for the general case of 
a Gaussian density matrix. In particular, for a purely chaotic system the two-particle cor- 
relation function must always be greater than one. On the other hand, in the presence of a 
coherent component the correlation function may take values below unity. The experimental 
situation is briefly discussed. 
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The most important method to obtain information about space-time aspects of particle 
production in high energy collisions is particle interferometry. One studies the correlation 
function of m (m > 2) identical particles, 

where k{ is the three-momentum of the i-th particle and Pi(k±, kg) are the ^-particle single 
inclusive distributions. For ultrarelativistic collisions, one usually considers bosons, i.e., one 
measures Bose-Einstein correlations (BEC). 

In practice most of the experimental measurements have been restricted so far to two 
particle correlations (see, however, ref. Jl|). Moreover, with very few exceptions, only 
identical charged particles were considered. Usually the data are analysed by fitting them 
to a form 

C 2 {kM) = 1 + \p(h-k 2 )\ 2 (2) 
where p(q) is the Fourier transform of the space-time distribution of the source elements, 

p(q) = J d 4 x p(x) exp(ig At x M ) (3) 

Thus, the inverse widths of the correlation function are related to the size and lifetime of 
the source. Eq. (0) was first derived by symmetrizing the two particle wave function and 
averaging over the emission points of the particles after taking the square of the matrix 
element 0. This corresponds to the assumption of a purely chaotic source (random phases). 
From Eq. (|2|) one obtains 

1 < C 2 (k h k 2 ) < 2 (4) 

and these are the bounds for C 2 in the conventional treatment of BEC. For neutral particles 
as well as for a partially coherent source, modifications of these bounds appear (cf. below). 
The existence of bounds for BEC is important, among other things, because they can serve 
as checks for results derived in any particular approach. 

A more general approach to BEC is based on the density matrix formalism, i.e., a formal- 
ism which does not assume pure states (wave functions). To be specific, we shall consider 
the case (which covers almost all phenomena known in quantum statistics) that the multi- 
particle system is described by a Gaussian density matrix []. If the particles are emitted from 
a large number of independent source elements, such a form of the density matrix follows 
from the central limit theorem (cf., e.g., ||). Upper bounds for the two-particle correla- 
tion function for a Gaussian density matrix have already been discussed in |EJ: for charged 
bosons, one finds again C 2 (ki,k 2 ) < 2, but for some neutral bosons like 7r°'s and photons, 
C 2 {k\, k 2 ) < 3. As will be demonstrated below, concerning the lower bounds one can derive 
the following general results: (a) for a purely chaotic source, C 2 (ki,k 2 ) > 1, and (b) for a 

1 We do not treat final state interactions here. 
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partially coherent source, (^(A/l,^) > 2/3 for identical charged bosons, while for 7r°'s and 
photons, C 2 (/ci,/c2) > 1/3. 

For a Gaussian density matrix, all multiparticle inclusive distributions can be expressed 
[§J in terms of the quantities 



D(k r ,k s )=JE r E s (a\k r )a(k s )) (5) 



D(k r ,k s )=-JE r E s (a(k r )a(k s )) (6) 



I(k r )=-i^/E r (a(k r )} (7) 

where a^(fc) and a(k) are the creation and annihilation operators of a particle of momentum 
k, and the indices r, s label the particles. 

For the general case of a partially coherent source, the single inclusive distribution can 
be expressed as the sum of a chaotic component and a coherent component, 

P x (k) = P1 hao (k) + P™ h {k) (8) 

with 

P° hao (k) = D{k,k) (9) 

and 

P?\k) = |/(fc)| 2 (10) 

The chaoticity parameter p, representing the ratio between the mean number of chaotically 
produced particles and the mean total multiplicity, is then 

D(k,k) 

P{k) = D(k,k) + \I(kW (U) 
To write down the correlation functions in a concise form, it is useful to introduce the 
normalized current correlators, 

[D(kr, K) ■ D(k a , k s )}2 [D(kr, kr) ■ D(k s , fc)]3 

where the indices r, s label the particles. Since d(k, k') and d(k, k') are in in general complex 
valued, one may prefer to express the correlation functions in terms of the magnitudes and 
the phases, 

T rs =T(k r ,k s ) = \d{k r: k s ) \ 
T rs =Tyk r ,k s ) = \d[k r , k s ) \ 

0rs=0 c/i (^V, k B ) = Arg d(k r ,k s ) (13) 
4> c ; h s =4> ch (k r ,k s ) = Arg d(k r ,k s ) 
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and the phase of the coherent component, 

= <p c (k T ) = ArgJ(Av). (14) 
The same notation will be employed for the chaoticity parameter, 

p r EE p(k r ) (15) 

For identical charged bosons (e.g., 7r~) the two-particle correlation function reads 

C 2 -(k u k 2 ) = 1+ 2 v /p 1 (l-p 1 )-p 2 (l-p 2 ) T 12 cos(0g-^ + ^) + piP2 3? 2 (16) 

For neutral bosons like photons, or 7r°'s, the terms d(k r , k s ) may contribute^ to the BEC 
function^: 

C™(h, k 2 ) = l + 2y/ Pl (l - Pl ) ■ p 2 (l - p 2 ) T 12 cos(^-^ + ^) + Pl p 2 T? 2 

+ 2yJ Pl (l-p 1 )-p 2 (l-p 2 ) f 12 C0S(#* - # - C 2 ) + p lP2 2* 

(17) 

Let us first consider the case of a purely chaotic source. Insertion of p(k) = 1 in Eqs. 
fliDD and ( pl) immediately yields C 2 (fci, fc 2 ) > 1. In the case of partial coherence, the terms 
containing cosines come into play and consequently C 2 may take values below unity. Eqs. 
( |H>| . |T7D imply that C 2 ~(k±, k 2 ) > 2/3 and C 2 °(ki, k 2 ) > 1/3. Because of the cosine functions 
in (|I6|J17D one would expect C 2 as a function of the momentum difference q to oscillate 
between values above and below 1. Indeed such a behaviour of the Bose-Einstein correlation 
function has been observed in high energy e + e~ collision experiments (cf., e.g., ref. 0), but 
apparently not in hadronic reactions. This observation was interpreted as a consequence 
of final state interactions in ref. |J. If final state interactions determine this effect, it is 
unclear why the effect is not seen in hadronic reactions. On the other hand, if coherence is 
responsible for it, this would be easier to understand. Indeed multiplicity distributions of 
secondaries in e + e~ reactions are much narrower (almost Poisson-like) than in pp reactions, 
which is consistent with the statement that hadronic reactions are more chaotic than e + e~ 
reactions 0. 

So far, two methods have been proposed for the detection of coherence in BEC: the 
intercept criterion [§] {C 2 {k,k) < 2) and the two exponent structure of C 2 0. Both these 
methods have their difficulties because of statistics problems and other competing effects. 
The observation of C 2 {ki, k 2 ) < 1 could constitute a third criterion for coherence, although 
C 2 (kx, k 2 ) < 1 is not a necessary condition. 

2 These terms are expected to play a significant role only for soft particles of energies of the order of the 
inverse lifetime of the source, or for sources with very small lifetimes B. 
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Recently, the two-particle correlation function has been calculated for photons emitted 
from a longitudinally expanding system of hot and dense hadronic matter created in ultra- 
relativistic nuclear collisions []IU |. For such a system, the particles are emitted from a large 
number of independent source elements (fluid elements), and consequently, one would expect 
the multiparticle final state to be described by a Gaussian density matrix. However, although 
the system is assumed to be purely chaotic (i.e., it does not contain a coherent component) 
the correlation function calculated in ITJ is found to take values significantly below unity. 
Clearly, this is in in contradiction with the general result derived above from quantum 
statistics (C 2 > 1 for a chaotic system). 

The expression for the two-particle inclusive distribution used in ref. Jl(| (equation (3) 
of that paper), which in our notation takes the form (w(x, k) is the emission rate) 



d 4 x 1 / d 4 x 2 w (x 1 , fci) w(x 2 ,k 2 ) [1 + cos((fci - k 2 )(x 1 - x 2 ))} (18) 



does not exclude values below unity for the two-particle correlation function^. Thus, the 
fact that the authors of ref. |lfj find values below one for C 2 may be due to the application 
of an inadequate expression for P 2 (ki,k 2 ) [] Q. 

— * — # 

Below, we shall briefly derive expressions for P 2 (ki,k 2 ) which could be used for the 



problem treated in |L0|, [IjJ and which do not suffer from the deficiencies mentioned above. 
To this end, we consider two conventional approaches used in this field: (i) the classical 
current formalism, and (ii) symmetrization of the wave function. 

(i) In the current formalism particle sources are described in terms of a distribution of 

HI 



classical currents l|T2l IT 



For a Gaussian distribution of currents and in the absence 



3 To see this, consider, e.g., the simple ansatz 

w(x,k) — const. exp[— a(x — /3k) 2 } S(t — to) 



where a and (3 are free parameters. The expression for P 2 (ki,k 2 ) used in rcf.jlTJ (cf. eq. fll8|)) then yields 



C 2 {ki,k 2 ) = 1 + exp 



q 

' 2a 



cos[Pf 



Clearly, if f3 exceeds a 1 the above expression will oscillate and take values below unity. On the other hand, 
in the current formalism (cf. below) one obtains with the same ansatz for w 



C 2 {k u k 2 ) 



exp 



9 

'2a 



> 1. 



4 Another possibility which was pointed out by the referee of the present paper is that the values below 
one obtained for C 2 in the calculations of ref. jIo| are caused by a certain approximation used to evaluate 
the integrals in JToj] . 

5 In a recent publication Jll[] the authors of ref. |10| have included the effects of transverse expansion. The 
resultant correlation functions do not take values below one. However, if the results for the 1-dimensional 
expansion arc affected by the use of an inadequate expression for P 2 {k\, k 2 ) 1 this would also cast doubts on 
the results obtained for the 3-dimensionally expanding system. 
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of a coherent component, all multiparticle distributions can be expressed in terms of the 
two-current correlator, (J*(x)J(x r )) [|J. The one- and two-particle inclusive distributions 
then take the formQ 

P x (£) = J <fx w(x,k) (19) 

and 

P 2 (k 1 ,k 2 ) = J d 4 xi J d^x 2 [w (xi, k\) w (x 2 , k 2 ) + w (x±, K) w (x 2 , K) exp [iq^x^ - x 2 )]] 

(20) 

where = (ki+k 2 )/2 and g M = k^—k 1 ^ are the mean momentum and momentum difference 
of the pair, and where 

w(x,k) = Jd A y (j*{ x + l)J{ x - V -)^ expi-ik^] (21) 

(ii) In the wave function approach, the amplitude for emission of two identical particles 
of momenta fci, k 2 from the space-time points x\, x 2 is 

A(ki, k 2 )=M[k l , x 1 )M(k 2 , x 2 ) exp[i{kix x + k 2 x 2 )} 

+M(k 1 ,x 2 )M(k 2 , Xl ) exp[i(A;iX2 + hxi)} (22) 
The one- and two-particle inclusive distributions are 

P^k) = J d*x \M{k,x)\ 2 (23) 

and 



P 2 (ki,k 2 )=- J d 4 xi J d A x 2 \A(ki,k 2 



d* Xl J d*x 2 [|M(fci,xi)r \M(k 2 ,x 2 )\ z 
+ M(k 1 ,x 1 )M*(k 2 ,x l )M(k 2 ,x 2 )M*(k 1 ,x 2 ) exp[i<^-^)]] (24) 
Both expressions [] ( |20| ) and ( p4"|) can be cast in the form 

P 2 (h,k 2 ) = P 1 (k 1 )P 1 (k 2 ) + D{k 1 ,k 2 )iy k {Kk 2 ) (25) 



6 For neutral particles, there may be an additional contribution to f^&ijfe) which only plays a role for 
soft particles and which will be neglected here, cf. footnote 2 on page 4. 

7 To obtain specific results for P2{k\,k2) for photons emitted from an expanding quark-gluon-plasma as 
discussed in jio), it is reasonable to substitute the emission rates |l5| for the integrands in Eqs. ([H]) and 
i-e., 

2.9 fc^f 



w(x, k) = \M(k, x)\ 2 = const. T 2 In ' JL" h 1 ) exp 



T 

where is the four velocity of the fluid, T the temperature and g the QCD coupling constant. 
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where for the current formalism 

h + k 2 



D(k 1 ,k 2 ) = J d^x w \ x, 1 - 2 j expl-ix^k^ - k%)] (26) 

and for the wave function approach 

D(k 1: k 2 ) = [ d 4 x M*(ki,x)M(k 2 ,x) exp[—iXfj,(ki — k 2 )\ (27) 



It then follows from Eq. (p5|) that C 2 (ki, k 2 ) > 1 both for the results of (i) and of (ii). 

The above considerations concerning bounds for the BEC functions refer to the case of 
a Gaussian density matrix. In general, a different form of the density matrix may yield 
correlation functions that are not constrained by the bounds derived here. For instance, for 



a system of squeezed states C 2 can take arbitrary positive values ||16|| . Moreover, for particles 
produced in high energy hadronic or nuclear collisions, the fluctuations of quantities such as 
impact parameter or inelasticity may introduce additional correlations which may destroy 
the Gaussian form of the density matrix and also affect the bounds of the BEC functions 



(cf. ref. ill, p 
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